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Abstract 

Wc show that both the mtcrior region r < M — \/ A'P — a? of a Kerr black hole and the a? > Af^ 
Kerr naked singularity admit unstable solutions of the Teukolsky equation for any value of the 
spin weight. For every harmonic number there is at least one axially symmetric mode that grows 
exponentially in time and decays properly in the radial directions. These can be used as Debye 
potentials to generate solutions for the scalar, Weyl spinor, Maxwell and linearized gravity field 
equations on these backgrounds, satisfying appropriate spatial boundary conditions and growing 
exponentially in time, as shown in detail for the Maxwell case. It is suggested that the existence 
of the unstable modes is related to the so called "time machine" region, where the axial Killing 
vector field is time-like, and the Teukolsky equation, restricted to axially symmetric fields, changes 
its character from hyperbolic to elliptic. 

PACS numbers: 04.50.+h,04.20.-q,04.70.-s, 04.30.-w 
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I. INTRODUCTION 

Kerr's solution [1] of the vacuum Einstein's equations in Boyer-Lindquist coordinates is 



where S = + cos^ 9 and A = — 2Mr + . Wc use the metric signature H to 

match the formulas in the original Newman-Penrose null tetrad formulation [2] and Teukol- 
sky perturbation treatment [3]. Kerr's metric has two integration constants: the mass M, 
and the angular momentum per unit mass a. They can be obtained as Komar integrals [4] 
using the time translation Killing vector field K"^ and the axial Killing vector field C,"' (in 
the above coordinates, these are d/dt and d/d(j) respectively). We will only consider the 
case M > 0, and we will take a > without loss of generality, since for a < we can always 
change coordinates — > —0, under which a — > —a. If < a < M (sub-extreme case), the 
E = ring curvature singularity at r = 0, ^ = 7r/2 is hidden behind the black hole inner and 
outer horizons located at the zeroes of A: = M — V — and Tq = M + V — a^. 
As is well known, and Tq are just coordinate singularities in (1), Kerr's space-time can be 
extended through these horizons and new regions isometric to I: r > Tq, II: rj < r < and 
III: r < r j arise ad infinitum and give rise to the Penrose diagram displayed in Figure 1. In 
the extreme case M = a, ri = To and region II is absent, however, we will still call region I 
(III) that for wich r > Vi = [r < Vi = To). In the "super-extreme" case a > M there is 
no horizon at all, the ring singularity being causally connected to future null infinity. This 
is not a black hole, but a naked singularity. 

In this paper we study the Kerr naked singularity, and region III of (sub-extreme and 
extreme) Kerr black holes. We will refer to these solutions of the vacuum Einstein equations 
as Kill and KNS from now on. These have a number of undesirable properties, among 
which we mention: i) the timelike curvature singularity as we approach the ring boundary 
at r = 0, 6 = it/2; ii) the fact that any two events can be connected with a future timelike 
curve (in particular, there are closed timelike curves through any point), making Kill and 
KNS "totally vicious sets" in the terminology of [5], and causing a number of puzzling 
causality problems [6]; iii) the violation of cosmic censorship. 

Two conjectures are known under the name of cosmic censorship: the weak cosmic cen- 
sorship conjecture establishes that the collapse of ordinary matter cannot generically lead 
to a naked singularity, while the strong cosmic censorship is the assertion that the maximal 
development of data given on a Cauchy surface cannot be generically continued in a smooth 
way [7]. The words in italics above signal aspects of the conjectures that need to be properly 
specified to turn them into a well defined statement. In any case, KNS violates weak cosmic 
censorship, and Kill violates strong cosmic censorship, since it is a smooth extension of 
the development of an initial data surface extending from spatial infinity of region I to spa- 
tial infinity of region I' in Figure 1, whose Cauchy horizon agrees with the inner horizon at r^. 

According to the Carter-Robinson theorem [8] and further results by Hawking and Wald, 
[9] if {Ai, gab) is an asymptotically-fiat stationary vacuum black hole that is non-singular on 




(r^ + a'^y - Aa^ sin'^ 6 
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FIG. 1. Penrose diagram for the maximal analytic extension of Kerr's space-time. Regions labeled 
I and r are isometric, and so are II and IF, and III , III' 

and outside an event horizon, then it must be an < member of the two-parameter 
Kerr family. The spacetime outside a black hole formed by gravitational collapse of a star 
is, independently of the characteristics of the collapsing body, modeled by region I of a Kerr 
black hole solution, placing this among the most important of the known exact solutions 
of Einstein's equations. Since this region is stationary (outside the ergosphere), there is a 
well defined notion of modal linear stability, under which it has been shown to be stable 
[3, 10]. It is our opinion, however, that proving the instability of those stationary solutions of 
Einstein's equations that display undesirable features, is as relevant as proving the stability 
of the physically interesting stationary solutions. In this line of thought, we are carrying 
out a program to analyze the linear stability of the most salient naked singularities (M < 
Schwarzschild spacetime, [11] Q'^ > Reissner-Nordstrom spacetime [12], and > 
Kerr spacetime [13, 14]; and also of those regions lying beyond the Cauchy horizons in 
Reissner-Nordstrom and Kerr black holes [12, 14]. As is well known, the linear perturbations 
of the spherically symmetric Einstein-Maxwell spacetimes are much easier to deal with than 
those of the axially symmetric Kerr spacetime, for which the only separable equations known 
to date are not directly related to the metric perturbation. They are equations satisfied by 
perturbations of the components of the Weyl tensor 

V'o := -CafeedrmW ^4 := -a^n^rnVm'^ (2) 

along a complex null tetrad i", n", m", m", among which the only non zero dot products are 
[2] 

rn„ = l, m"m, = -1. (3) 

We use a bar for complex conjugation, l"" and are real vector fields, whereas is complex. 
The null tetrad we use is that introduced by Kinnersley [15], given in equation (4.4) in [3]. 
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If we take 

Sabcd = ? 4! l[ani,m,rhd] (4) 
as a right handed volume element, we find that the following two-forms are self dual 

■m[anb], n[ah] + 'm[a'fhb], /[^mb] (5) 

A complex electromagnetic field can be written as 

Fab 20i(n[aZ6] + m^afrib]) + 202 ^[a"^6] + 200 fh[anb] 

+ 201 {n[alb] + rfi^arrib] ) + 202 l[afhb] + 20o m[„nb] . (6) 

If Fab is real then 0^ = 0^, if Fah is self dual (anti-self-dual) then the 0j (0j) vanish. Teukolsky 
[3] found that Maxwell, (Weyl) spinor and scalar fields on a Kerr background can be treated 
in a similar way as the gravitational perturbations ^■^o and 5'04. If we take the components 
of the Maxwell fields (see (3)) 



Fabl^'m' 



Fabitn' + m'^m') 02 = F„6mV^ 



(7) 



and those of the two component spinors xa 



Xo = Xao 



Xl = XAi , 



(8) 



and weight them with an appropriate power of the spin coefficient p — m^wf'Vbla — 
{iacos9 — r)~^: 

*i := Xo, := P~\u *i 0o, P'^h: *2 := (^V'o, *-2 := P'^Si/j^, (9) 

then the (source-free) Maxwell, spinor, and linerized gravity equations can all be encoded 
in Teukolsky's master equation [3] 



(r^ + a2)2 



A- 1- ( A 

or 



dr 
-2s 



1 d 



+ 



sin 9 89 

M{r^ - a2) 



A dtdcj) 



sin 9 



09 

— r — ia cos 9 



-2s 



A sin^^J 902 
a(r — M) i cos i 



sm 



d(j) 



dt 



+ (s2cot2^-s)*, = 0. (10) 



The index s in gives the spin weight under tetrad rotations. The above equation also 
gives the massless scalar field equation D^o = if we set s = 0. 

In [13] we found numerical evidence that there are solutions of the s = —2 Teukolsky 
equation in the KNS that grow exponentially in time while satisfying appropriate boundary 
conditions. In [14] we confirmed this fact by proving that there are infinitely many axially 
symmetric unstable (meaning, behaving as e'^* for some positive k) solutions of the s = — 2 
equation in the KNS, and also in Kill. 

In this paper we extend further this result to other linear fields. We show that there are in- 
finitely many unstable solutions of the Teukolsky equation for any s value (|s| — 0, 1/2, 1, 2). 
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The proof is given in Section III, with some calculations relegated to the Appendix. The 
existence of the unstable modes is shown in Section IV to be related to the time machine 
region near the ring singularity, which produces a change of character of the Teukolsky PDE 
-restricted to axial modes- from hyperbolic to elliptic. It is suggested that the emergency 
of an instability when a PDE changes from hyperbolic to elliptic is generic, and this is 
illustrated with a simple toy model in 1 -|- 1 dimensions. 

The use of unstable solutions of Teukolsky's equation as "Debye" potentials [16] for con- 
structing unstable spinor. Maxwell and linear gravitational fields is illustrated in Section V, 
where the reconstruction process is explained in detail for Maxwell fields. All the recon- 
structed fields decay properly along spatial directions whilst growing exponentially in time. 
This is the notion of instability used in this work. The lack of a sensible initial value formu- 
lation due to the fact that there are no partial Cauchy surfaces in Kill and KNS, forbids a 
more traditional approach to the stability issue. This is quite different from what happens 
for the Reissner-Nordstrom and negative mass Schwarzschild timelike naked singularities, 
for which a unique evolution of data given on a partial Cauchy surface can be defined, and 
instability proven afterwards [11, 12]. 

The following section contains information on the Teukolsky equation that is used in the 
proof of existence of unstable modes. 



II. SEPARATED TEUKOLSKY EQUATIONS 

Introducing 

= Ro,,m,s{r)SZ',{9) exjp{tmcl)) cxp{-iut), (11) 
the linearized Teukolsky PDE (10) is reduced to a coupled system for S and R, 



1 d f.^^ndS\ , / 9 .2 9z) o„, {m + scos9) 



2 



sin^— + aWcos^e- 2aujscos9 - — + E - ) S ^0 (12) 



sine d9 V d9 \ sin^ 9 



^(fR ^ ^^dA dR 

dr^ dr dr 

( - 2is(r - M)K ^ ^ 22/ m1 ^ ^ / ox 

+ < ^ ^ h 4:irujs -[E- 2amu + a^u^ - s{s + 1)] ^ = 0, (13) 

where K = [r"^ + a'^)u — am. The system (12)-(13) is coupled by their common eigenvalue 
E, whose relation with the separation constant A in [3, 17] is given hj A = E — s{s + 1). 
Suppose s, m and cu are given, then E in (12) has to be chosen so that S is regular on the 
sphere. This gives a denumerable set of eigenvalues that we label as £'""'^(s, m, acu), with 
i = 0, 1,2, ... and E increasing with £. In a similar way, E in (13) is chosen such that R 
decays properly as |r| 00 for the KNS (r — )■ —00 and r — )■ r~ for Kill), and this also 
gives a denumerable set of increasing values E^^'^{s,m,aio),n = 0,1,2,.... A solution of 
the system (12)-(13) is obtained whenever E^"^{s,m,au!) — E'l^'^{s,m,au) —: E. Thus, for 
given (s, m), we may regard a solution as an intersection of the curves vs au and 
vs acj, the allowed frequencies being those at which the curves intersect. This point of view 
is the one used in the proof of instability below, for which we restrict our search to auj = ik, 
k a real positive number (so that (11) gives an exp{kt/a) behaviour), and show that there 
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are intersections for every i and, at least, the fundamental radial mode n = 0. 

One way of finding the radial and angular spectra consists in reducing the regularity 
conditions to a continued fraction equation involving E [18]. This equation arises when 
solving a three term recursion relation on the coefficients of a series solution for S and K in 
(12) and (13) [18]. An alternative way to obtain the angular spectrum, which is well suited 
to the case we are interested, m = and aio = ik, is discussed in Section 2.1 of [14]. This 
approach is used in the numerical computations leading to Figure 2 below. 

A. Angular equation: spin weighted spheroidal harmonics 

The solutions of (12) which are regular on the sphere arc called spin weighted spheroidal 
harmonics (SWSH). The spectrum of E values is discrete, and we will use the notation 
E^"'^{s,m,au)), i = 0, 1,2, ... to label the eigenvalues in increasing order for a given set of 
parameters {s, m, acu) (the notation is not unified in the literature, note that we use £ — 
for the lowest eingenvalue independently of the spin weight s). Equation (12) exhibits some 
interesting symmetries: if {S{6),E) is a regular solution of (12) for some given {s,m,u) 
values, it is easy to check that (5'(7r — 9),E) is a solution for (s, — m, —cu) and for (— s, m, cu), 
and that {S{9),E) is a solution for {s,m,u). This implies that 

E'^"'^{s, m, acu) = E'^"'^{s, —m, —aoj) = Ef^^{—s, m, aou) = E'P^i^s, m, aCj). (14) 

We are interested in axially symmetric (m = 0) solutions with purely imaginary frequen- 
cies {aoo = ik,k & M.). In this case, from a solution S{6) we can extract solutions with real 
and imaginary parts of opposite parities by taking the linear combinations S{9) ± S{n — 0). 
Also, the eigenvalues are real, as follows from (14): 

E'P'is, m = 0, ik) = EP{s,m = 0,tk). (15) 

The behaviour of E'^"'^{s,m = 0,auj = ik) in the limits A; ^ 0"*" and k ^ oo will be rel- 
evant in what follows, since, as explained above, an intersection of the curve E'^"'^{s,m — 
0,auj — ik),k e M"*" with any of the radial curves implies the existence of an unstable 
(i.e., ~ exp{kt/a)) mode. To study these limits we will assume that s > without loss of 
generality (see (14)). Setting x := cos{6),m = 0,auj = ik in (12) we find that, for any k, 
this equation has regular singular points at a; = ±1, the possible behaviour of local solutions 
around these points being (1 — |a;|)*/^ or (1 — |a;|)~'^/^. 

For A; = 0, we expand the regular solutions as 

oo 

S{x) = (1 - a;')(^/2) - ^)'' (16) 

and find that (12) implies the recursion relation 

2(j + l)(j + 1 + s)a,+i - ((j + s)(j + s + 1) - E) a, = 0, (17) 

which, for large j, gives a^+i ~ %/2. The series in (16) will thus diverge at a; = — 1 unless 
we cut it to down to a polynomial of degree ^ = 0, 1, .. by choosing E = {i + s){i + s + l) for 
some £ — 0,1, 2, .... Repeating the calculation for negative s, or just using (14), we obtain 

E ^ E^'^%s,m^O,auj^O) = {£ + \s\){£ + \s\ + 1). (18) 
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k k 

FIG. 2. The left panel shows £"""^(5 = 1/2, m = 0, aw = ik) obtained numerically by solving eq. 
(9) in [14] for s = 1/2 and different values of G [0,15], together with a least square linear fit 
using the large k data points, which gives E = 2.99843fc+constant, in excellent agreement with the 
expected asymptotic expansion Eq. (20) . The right panel shows the values obtained in this way for 
< k < 2.5. The solid line is the low frequency approximation of £'^"^(s = 1/2, m = 0,auj = ik) 
in [17, 19] to order k^. 

A more detailed analysis of (12) using continued fraction techniques gives a Taylor expansion 
for Ef^^ {s,m,au) for complex uj near cj = The expansion up to order (ao;)^ is available 
in the literature (see ([17, 19] and references therein), exhibits the symmetries (14) and has 
(18) as the leading order term, i.e., 

E = Ef'\s,m = 0,acc;) = (£ + |s|)(£ + |s| + l) + 0{u). (19) 

Asymptotic expansions for au = ik,k ^ 00 can be found in [17, 20-22]. Particularly 
useful to our purposes is that, in our notation [17], 

Ep%s,m,au ^ik) ^ (2e + l)k + 0(k°), as A; ^ 00 (20) 

We should warn the reader, however, that a complete proof of the above formula is not 
available for s 7^ 0. Although arguments suggesting the vahdity of (20) are given in [17], 
where the formula was also numerically checked for s = 1, 2 the case s = 1/2 has not been 
reported as tested there. Given that these equations are key in the proof of instability for 
spinor fields, we tested numerically their validity using the method developed in Section 2.1 
of [14]. We have found and excellent agreement with both (18) and (20) for s = 0, 1/2, 1, 2. 
As an illustration, we give in Figure 2 the results for the only case not dealt with in the 
literature, that of s = 1/2. 

B. Radial equation: reduction to a Schrodinger form 

Equation (13) is of the form AR + QR + {Z — E)R = 0, dots denoting derivatives with 
respect to r. If we introduce an integrating factor L, ip :— R/L, and change the radial 
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variable to r*, where ^ 



: J, with / an unspecified positive definite function of r, (13) 
gives the following equation for ip: 



Lf 



L" 

T 



QfL' fz 

LA A 



f2 



(21) 



where primes denote derivatives with respect to r*. By choosing / = VA (note that A 
is strictly positive for Kill and KNS) and L such that the coefficient of ip' vanishes, i.e., 



.(Hf+li 



L — A 4 , (21) reduces to a stationary Schrodinger equation with energy eigenvalue —E, 



Ui/j := -i/j" + Vij^ -Ell) 
In the case m — 0,auj — ik, the potential is 

"r (r^ + ra^ + 2Ma^ " 



(22) 



a2 (r2 -2Mr + a2) 



r+2s 



.3 _^r2j^ j^^^2 j^^2 



+ 



1 + 



(M^-a^)(4g^-l) 
r2 - 2 Mr + a2 



a (r2 - 2 Mr + a^) 
A;V2 + A;^! + (23) 



k 



rad I 



We will show that the spectrum of (22) is entirely discrete, and use the notation —E^^ 
0, auj — ik),n — 0,1,2, ... for its eigenvalues, to be consistent with our previous conventions. 
For k large enough, V is negative in a region r„ < r < 0, and the resulting bound states 
lead to unstable modes of the Teukolsky PDE. This is explained in detail in the following 
Section. 



III. UNSTABLE LINEAR FIELDS ON KERR' S SPACETIME 

The results from the previous Section can be summarized as follows: there are solu- 
tions of the Teukolsky equations behaving as e*^*/" if and only if solutions can be found of 
equations (12) and (13) with the same E value for u = ik/a, i.e. E^^^{s,m,LO = ik/a) = 
£^™°'(s, m,uj = ik/a) for some £ and n. Since we restrict our attention to the axially symmet- 
ric case, we will drop the m index from now on, and use (22)-(23) instead of (13). Given that 
the instability is a consequence of the intersection of spectral lines of the angular and radial 
operators for purely imaginary frequences, we need to gather information on the spectra of 
theses operators. Since we are interested in spotting intersections for cu = ik/a, k & (0, oo), 
we will gather information on the asymptotic expressions for these spectra in the limits 
where A; — >■ 0+ and /c — )■ oo. The strategy of the proof consists in showing that in one of 
these limits E'/^^ > E"™*^ whereas in the other E^""^ < Eq°''^, thus the intersection follows 
from continuity on the spectral lines on k. For the angular equation, these limits are given 
in (19) and (20). For the radial equation, we need to work them out and, since the analysis 
depends on the domain of r and boundary conditions, we will consider separately Kill and 
the KNS. 



A. Unstable modes on a Kerr naked singularity 

In this Section we consider the extreme case > M^, for which — oo < r < oo,t,9 and 
(f) are global coordinates, and A > everywhere. The choice ^ '■= j = I/'n/A made above 
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gives an adimensional r*, 



, ,r-M+Vr^-2Mr + a^\ fin r ^ oo 
M r^-co 



that grows monotonically with r, and can easily be inverted in terms of elementary functions, 

'^^ + M + — -— . (25 

Note that (22) defines a quantum mechanical problem in the entire r* line, with H in (21) a 
self adjoint operator in L^(M, dr*). The asymptotic form of the potential (23) for large \r*\ 
does not depend on the value of s, and is given by 

V-^< fiM^-a^)kV _2.* . . (26) 



2Ma 



e r ^ — oo. 



Prom the above equation and the fact that V is smooth we conclude that V reaches a mini- 
mum and H in (21) has an entirely discrete, bounded from below spectrum —E'!^'^{s, k),n — 
0, 1, 2... As explained above, we will need information on this spectrum both in the A; — > 0"*" 
and k ^ oo limits. 



1. Asymptotic behaviour of the radial equation spectrum as k oo 

The large real k limit is simple to deal with, because the behaviour of V in this limit 
does not depend on the value of s, and, therefore, the analysis in [14] for s = —2 applies 
with only minor modifications. The cubic polynomial + a^r + 2Ma^ in the numerator of 
V2 (see (23)) has a unique real root at r = r„(M) < 0, thus V2 is negative in the interval 
r„(M) < r < 0, and nonnegative elsewhere. Note that since = — r„^/(2M + r„), then 
rn{M) goes from —M to —2M as goes from to infinity. 

Let ip he a. properly normalized — dr* — 1) function supported in the 

interval (r„(M),0), then 

2 

{mm = (V'l - (9/9r*)» + kHmm (27) 

j=0 

with 

dr 

{mm = m'V,^<0 (28) 

Take kc to be the largest among zero and the real roots (if any) of 

p{k) - (V'l - (9/ar*)» + (v^i^oiV') + mvim + y(V'|^2|V') 
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(note the one half factor in the k'^ term!), then p{k) < ior k > kc and, if —Eq^'^{s, aui — ik) 
is the lowest eigenvalue of H (see equation (22)), 

-Eo™'(5,a^ = tk) < imi^) = y(V^I^2|V^) +P(^) < y(V^|^2|V') if ^ > k,. (29) 

meaning that the absolute value of the fundamental radial level grows at least quadratically 
in A; ^ 

E'^-^{s,aw = ik) > ^\{i;\V2m , if A; > A;e (30) 



2. Asymptotic behaviour of the radial equation spectrum as k 0'^ 

A quick inspection to the potential (23) gives the the minima for A; = and different spin 
weights: 



i ,s = 



mm{V{r, A; = 0, s), r G M} = < 
A few lengthy calculations show however, that 
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i ,1^1 = 1 

|s| = 2 



_7 
2 ' 



lim min{y(r. A;, s, a > M), r e R} = 7 - s^ (32) 

fc— >-o+ 4 

so that there is a discontinuity at A; = for higher spin values. Since we will be using 
continuity arguments in our proof of instability, wc will consider the fundamental energy 
—Eq"''^{s, au! = ik) of the quantum Hamiltonian in (22)-(23) as a function of A; e (0, 00), for 
which 

lim -El'"^{s, aw = ik) > lim min{y(r. A;, s), r e R} (33) 

fc->-0+ fe-)-0+ 

given in (32). Appendix A contain the details of the calculations leading to (32). 
3. Proof of the existence of unstable modes for every s 

Let us gather the relevant results of the previous sections for the axially symmetric 
(m = 0) modes. For A; ^ 0+ and £ = 0, 1, 2... 

£;p(s, aw = ik) \k=o+= {i + \s\){e + |s| + 1) > - i > E'o"\s, au = ik) \k=o+, (34) 
whereas for large enough positive real k 

E^"^{s,au;^ik)^2{e + l)k + 0{k'^) < —MV2\^)\ < El'"^{s,auj ^ ik). (35) 

By continuity, we must have, for every I and s, a A;(£^s) such that E'^"'^{s,au = ik(^e,s)) = 
£'^'^(5,00; = ik(^e^s))- This proves that there is an axially symmetric unstable solution of 
the Teukolsky equation for the fundamental radial level and every harmonic number i. For 
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higher excited radial level the arguments in [14] for s — —2 suggesting that there also are 
intersections generalize to arbitrary s. In any case, we have shown that there are infinitely 
many unstable modes for every spin weight. These solutions of the Teukolsky equations 
decay exponentially with |r| as |r| — )■ oo, so that they are initially bounded, and grow 
exponentially in time. 

The calculations above can be adapted to perturbations in the interior region r < ri :— 
M — V — of an a < M Kerr hlack hole. There are some subtle differences between the 
extreme a — M and sub-extreme a < M cases, as shown in the following sections. 



B. Unstable modes on region III of an extreme Kerr black hole 



For the extreme black hole the solution of dr*/dr = l/v A in the interior region r < r. 
To = M is 

'M - r 



r* = -ln(^^^J,r<r„ (36) 
with inverse 

r = M(l -e"''*), -oo<r*<oo (37) 

Using the integration factor A ~ as before, we are led back to (22) and (23), with r given 
in (37). Note that 

^ ^ f4/c2exp(2r*) ,r*^oo 

|fc2exp(-2r*) ,r*^-oo, ^ ' 

then for any A; > the spectrum of the self-adjoint operator "H is again fully discrete and 
has a lower bound. The argument leading to (30) goes through in the super extreme case 
without modifications, because the test function in (27) is supported in the r < region. 
Thus, the fundamental energy of the radial Hamiltonian is negative and there is an such 

that El°''^{s,au = ik)\k=o+ < {io + |s|)(^o + \s\ + 1), from where it follows that there is an 
unstable mode for every £ > The radial decay of these modes as r ^ and r — oo is 

^^|(llfc)"""*exp[-2i(Ji,)](H-0(MF-)) ,r^M- 
\{f exp [g] (1 + 0(M/r)) ,r ^ -oo 

C. Unstable modes on region III of a sub-extreme Kerr black hole 

In the sub-extreme case dr*/dr — 1/VA. and r < give 

^i + ^o-2r-2v/(ro-r)(r,-r) . 

' n + To ' 



so that r* has an upper bound: 



ro + ri 



-oo < r* < r* In ( ^^^^ ) . (41) 
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Near the domain boundaries, 



(42) 



and 

|[A;(M2 -a2)/(2aM)]2 exp(-2r*) ,r*^-oo ^ 

where we have defined 

/ T • / 'T -I- T" - \ 

- s (44) 




The sub-extreme case is essentially different from the extreme and super-extreme cases 
because (22) is a Schrodinger equation on the half axis r* < r*. with a potential that is 
singular at the r* boundary. This situation and type of singularity is well known [23, 24]. 
Any local solution of (22), in particular those which are square integrable for r* near — oo, 
behave as 

V' ~ a Ur* - r*)^+" + ...1 + b Ur* - r*)^"" + ...1 , (45) 

near the horizon. Thus, if > 1, these are not square integrable near the horizon, unless 
6 = 0, and this is precisely the condition that selects a discrete set of possible E values as 
the spectrum of T-L, and that defines de space of functions where TL is self-adjoint. This case 
is called limit point in [24]. It is quite different from the limit circle case z/ < 1, for which 
for any E the eigenfunction behaving properly at minus infinity will be square integrable 
in r* e (— oo,r*), and a choice of boundary condition needs to be imposed to define a set 
of allowed perturbations Dphyg, in order that ?^ be a self-adjoint operator on Dp^ys (i.e., 
Dphys = D*f^y^, see [24] for more details), and thus have a complete set of eigenfunctions. 
This is done by requiring a behaviour like (45) with a fixed (possibly infinite) b/a ratio [24]. 
Since we are ultimately interested in the large k case, in view of (44), we do not have to 
deal with this ambiguity. In any case, regardless of our choice of boundary conditions, the 
test function used in (27)- (28), being supported in a r < region, will belong to the chosen 
space of perturbations, and the argument of instability used for the nakedly singular Kerr 
spacetime will go through in the sub-extreme black hole case if, as done for the extreme 
case, we restrict the harmonics to i > ig with io the smallest non-negative integer satisfying 
£;7''(s,aa; = iA;)|fc=o+ < (4 + k|)(4 + |s| + 1). 



D. Consistency with modal stability outside the black hole horizon 

We would like to note that our results do not contradict the well established modal 
stability of the outer stationary, region I of Kerr black holes (see [10] en references therein). 
Our arguments break down in this case since the interval where V2 is negative hes outside 
the domain of interest. 

Consider the extreme case a — M, and switch to the adimensional variable x — r/M, then 



Vext — 



x{x + l)(x^ — X + 2) 



[X 



k^ + 2s 



X 



2x-l 



X 



(46) 
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For s = 0, Vext > 1/4 outside the horizon {x > 1) and thus there is no instabihty. For s 7^ 0, 
dVext/dx = at Xo gives 

k {xq — ^){xo — + 3) fA7\ 

~s ~ {xo'' + l){xo^-2xo-iy ^ ^ 

The function on the r.h.s. above decreases monotonically from zero to minus infinity for 
Xo e (1,1 + V2), then from infinity to zero if x e (1 + -\/2, 00). Thus, for any s 7^ 0, Vext 
has a unique critical point, which (by inserting (47) in d'^Vext/dx'^) we find that is a local, 
then absolute (using that Vext 00 for a; — )■ l"*" and x — )■ 00) minimum V^^f- of Vext in the 
domain of interest. This absolute minimum V°^^ can easily be seen to decrease with k, with 
a lower bound | — as A; — > 0"^. Summarizing 

Vext{k,x,s)>^- s'',\s\^ 0,1/2,1,2; k > 0, x > 1. (48) 
Then for the radial equation we get 

El'"^{s,aw^ik) < - ^ < + |s| < (£+ |s|)(£+ |s| + 1) < Ef'%s,au ^ ik), k>0,l>0 

and thus the intersection argument implying the existence of modes that grow exponentially 
in time breaks down. 

The reasoning for the sub-cxtrcmc case is similar, although the calculations are more com- 
plicated. Instead of (47) we get a two branched solution for k, only one of which corresponds 
to local minima. The bounds (48) are obtained again, and thus there is no instability. 



IV. INSTABILITIES AND TIME MACHINE 

The axial KiUing vector field of Kerr spacetime (C = d/dcf) in Boyer-Lindsquit coordi- 
nates) becomes timelike in what is called the time machine region ._5^of Kerr spacetime [5] , 
the region where 

(r^ + cos^ 6) (r^ + a^) + 2Mra^ sin^ ^ < 0. (49) 

The character of the restriction of the Teukolsky PDE (10) to the space of functions satisfying 
= im^) changes from hyperbolic to eUiptic within To see this, note that the second 
order terms of Teukolsky equation are independent of s, and so, for any s value, they equal 
those of the scalar wave equation 

= 5"'V„V,^ = -^da {yW\g'^'d,^) ^ g'^'dad,^ ^ (^g""' - 8^8,^ (50) 

where ~ means "equal up to lower order terms" , and £(^^ — im^ was used in the last step. 
The proof of existence of unstable modes in Teukolsky equation in the previous Section is 
based on the fact that the piece V2 of the potential, which is dominant for large k, becomes 
negative in the region r (r^ + ra^ + 2 Ma^). This region is precisely the intersection of ^ 
with the equatorial plane = ti/2 (see (49)). 

The existence of instabilities seems to be related to this change of character of the Teukolsky 
PDE for axial modes from hyperbolic to elliptic. To illustrate this point, consider the simple 
toy model {a,uio positive): 

ujHx' - a') 8^^ 8'^ 
4 9^-9^=°' 
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which is hyperbohc for \x\ > a, eUiptic otherwise. Unstable solutions ^{t,x) — e'^*/"'0(x) 
will exist if 

'^^ + f^y(,2_„2)^^0. (52) 



The above equation is that of a quantum harmonic oscillator, it has square integrable solu- 
tions if 

E = kuo (n+l^- ' = 0, n = 0, 1, 2, .... (53) 

Thus, there are no instabilities if a = (i.e., (51) is hyperbolic everywhere). Otherwise there 
will be infinitely many unstable modes, with 

One of the consequences of the existence of the time machine region is that it allows to 
construct a future directed timelike curve connecting any ordered pair of events in either 
Kill or KNS (sec [5] and references therein). This causes a number of difficulties when 
trying to define notions such as "evolution" and "instability" , as discussed in the following 
Section. 



V. UNSTABLE MODES AS DEBYE POTENTIALS 

Solutions of the transposed Teukolsky equations of different spin weights can be used as 
"Debye potentials" to generate Maxwell, spinor, and linearized gravity fields [16]. An expla- 
nation of why this is so was first given by Wald in [25] , and is reviewed in detail in appendix 
C of reference [26]. It is based on a notion of transpose of a linear differential operator O 
acting on tensor fields of rank /c, under the inner product (t/, V) := Jj^ t^ai...afc^"^ '"*, where 
M is the spacetime manifold and indices are raised and lowered using the metric. The 
transpose is defined as usual by {U,OV) — {0'^U,V), and assumes a proper decay of the 
fields in the domain of O, so that integration by parts is allowed. Suppose / is the tensor 
field we are interested in (e.g., the Maxwell potential Ah, or the metric perturbation hab), 
S{f) = the linear differential equation that it satisfies. In the Teukolsky formalism, with 
the exception of the scalar field equation, one does not work with the field / of interest, but 
with a derived field ips '■= %{f)- Here 7^ is a linear differential operator that projects out a 
null tetrad component with spin weight s of a tensor derived from /, e.g, a perturbed Weyl 
component associated with a metric perturbation hab- Teukolsky perturbative treatment 
can be summarized as follows [25]: there exist hnear differential operators Sg and Og such 
that 

S,£if)^0,rs{f)^OsM. (54) 

The field equation E{f) — then implies Teukolsky equation OgO's) — 0. The operators Og 
are defined by the left hand sides of the following equations in [3] (refer also to equations 
(2)-(8) above): (2.12) for s = 2 {^2 = Sipo), (2.14) for s = -2 (^_2 = S1/J4), (3.5) for 
s = l {iPi = 0o), (3.7) for s = -1 (^_i = 02), (B.4) for s = 1/2 {^y2 = Xo) and (B.5) 
for s = —1/2 {ijj-i/2 = Xi)- Using the information in Table I of [3] (and calculating for the 
spinor case), we find that the relation between the Os and the operator Tg in the master 
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Teukolsky equation (10) are 



f(2S)-oT. ,.>0 
' \(2E)-V-'^oT,op2. ^s<0 ^ ^ 

Thus, Os{ips) = reduces to Tg'^g = 0, where ^1/^ = ipg for s > 0, and = p^'''?/'s for s < 
(cf. equation (9)). Teukolsky master equation T^^^ = is spelled out in Boyer-Lindquist 
coordinates in (10). 

As Wald noted in [25], for spinor. Maxwell and linear gravitational fields, — E, the 

transpose of iS^^ = OgTs (equation (54)) then gives %^Os^ = £Ss^ ■ Thus, if ipg is a solution 
of the transposed Teukolsky equation, Og'^ips = 0, then SS.J'ips = 0. In other words, ips is 
a "potential" for a solution / = Ss'^ips of the field equation £{f) = 0. 

A straightforward computation shows that there is a close relation between transpose and 
spin weight flip: 

p^NoO(±/op-2N = 0(^,) (56) 

Solutions for the transpose of the spin weight s source free Teukolsky equations can then 
be readily obtained by multiplying a solution of spin weight —s times an appropriate power 
of p. Unstable solutions of the Teukolsky equations will therefore produce unstable spinor. 
Maxwell or gravitational flelds, since the exp{kt/a) factors in the potential go through the 
differential operators Sg^. 

When analyzing the linear stability of a super-extreme Reissner-Nordstrom spacctime (or 
the interior static region of a Reissner-Nordstrom black hole) one is faced with the problem 
that the unperturbed spacetime is non globally hyperbolic due to the timelike singular 
boundary [12]. The evolution of fields on this spacetime is a priori not well defined, and the 
curvature singularity poses the additional problem of deciding what should be considered 
a "reasonable" behaviour for linearized perturbations. The way around these problems 
is hinted by the observation that there is a unique choice of boundary condition at the 
singularity that guarantees that the perturbed curvature invariants will not diverge faster 
than the unperturbed ones as the singularity is approached [12]. By choosing this partic- 
ular boundary condition we make sure that the perturbation treatment is self- consistent, 
as perturbations can be uniformly bounded on an "initial time" partial Cauchy surface 
Eo that meets the singularity (any hypersurface orthogonal to the timelike Killing vector 
field). At the same time we solve the issue of uniqueness of evolution from data given at Sq. 
Furthermore, this evolution preserves the chosen boundary condition [12]. KNS, as well as 
KII, also have a timelike curvature singularity, the ring singularity, located at r = 0, 6' = 7r/2 
in Boyer-Lindquist coordinates. Note, however, that r e (— oo, oo) for the super-extreme 
case (r G (— cxD,rj) for the black hole interior), as one can enter the r < region avoiding 
the singularity. The character of the singularity and the chosen fields is such that the 
separated Teukolsky equations (12) and (13) are not singular, i.e., r = is a regular point 
of the ODE (13), and similarly for 9 — 7r/2 in (12). When solving (13), which is a second 
order equation, one can impose that R{r) vanish as r — > — oo and r — > oo (r — > rj), but 
that leaves out any further choice, such as a selecting a specific behaviour as r — )■ 0. This 
implies that, although (super extreme or black hole interior) Reissner-Nordstrom and Kerr 
spacetimes share some properties, such as the lack of a Cauchy surface and the existence of a 
time-like singularity, the issue of field propagation on those spacetimes is technically rather 
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different, the Kerr ring singularity being milder. On the other hand, the causality issues are 
much worse in the Kerr case. This is because, as mentioned in the Introduction, any two 
events in KII or KNS can be connected with a future directed timelike curve (in particular, 
there is a closed timelike curve through any point.) There is no partial Cauchy surfaces 
and thus no clear notion of "initial time slice" that allows to pose the stability problem as 
an initial value problem. The t = constant slices are spacelike outside a compact set, and 
our notion of instability is limited to the observation that there exist solutions to the linear 
field equations behaving as exp{kt/a), k > 0, decaying exponentially or faster as |r| — >■ oo in 
KNS (vanishing at the inner horizon in Kill), and behaving "properly" -as defined below- 
near the ring singularity. 

To check how unstable fields behave near the ring singularity of Kerr spacetime we may 
use unstable solutions of the Teukolsky master equation (proved to exist for every s in 
the previous Section) as Debye potentials for unstable spinor. Maxwell or gravitational 
fields. Note, however, that, in contrast to our previous instability results for the Reissner- 
Nordstrom or naked Schwarzschild case case [11, 12], we lack of explicit expressions for the 
unstable solutions of the Teukolsky master equations. We can still get information on the 
behaviour of unstable fields near the ring singularity by analyzing the Probenius series so- 
lutions of the ODE (13) near r = and the ODE (12) near 9 = tt/2. These series, however, 
are independent of the stable or unstable character of the solution, since u does not show 
up at leading order. Thus, whatever criterion we adopt to disregard field solutions from 
Debye potentials based on their behaviour near the singularity, it will overrule every field 
(unstable or not) that can be constructed using the potential method outlined above. 



A. Maxwell fields 

For Maxwell fields s = ±1, and the operators and fields in (54) are 

/ = A,, [E{A^)\a = V^VeA„ - V^V„A„ (57) 

and, for s = 1, 

TMb) = l''m\VaA, - VkAa) (58) 
S-,{Ja) = \{5-l3-a-2T + 7:){jr) -\{D~2p- p)U,m') (59) 
Oi{iPi) ^{D-2p-p){A + ii- 27)V^i -{5- /3-a-2r + 7r){S + 7r- 2a)V^i, (60) 

where the standard null tetrad formulation notation [2, 3] is used {D, A, S, 5 are derivatives 
along the tetrad vectors, the other symbols represent spin coefficients.) Since 7i projects out 
a self-dual piece of Fab (see (58)), the complex potential Si^ipi constructed from a solution 
Ci^(V'i) = of the transpose Teukolsky equation [25], 

[Si'^^i]b = [-h{S + 2P + r)+ mb{D + p)]Vii, (61) 

will produce a self dual Maxwell field Gab = Fab + i*Fab [25]. In fact, the easiest way to 
check that the exterior derivative Gab of the potential (61) satisfies the source- free Maxwell 
equations is by checking that it is self-dual, which amounts to checking that the contractions 
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of Gab with any of the three anti-self-dual two-forms obtained by complex conjugation of (5) 
vanishes as a consequence of Oi'^{iIji) = 0. 

To evaluate the strength of the real Maxwell field Fab near the ring singularity, we compute 
the algebraic invariants h — FabF"-^, I2 = Fab*F°'^ (any other algebraic invariant will be a 
polynomial on these). Note that, since ^GabG"-^ — Ii -\- i/2 I, we can compute the 
invariants of Fab niore efficiently without even taking the real part of Gab- For generic 
separable solutions = e^^^R{r)S{d) of the s = —1 Teukolsky mater equation, (61) gives 
a field whose invariants admit an expression that can be simplified near the ring singularity 
by applying iteratively the equation T^ilipi) = 0, to 



2 f-iS{e)a^ + R{r)^) e2-*/a 

^ ^ ^ r^- — 7^ ■ (^2) 

(r -I- la cos [B)) 

As already explained, this leading order term (omiting the exp{2iu)t/a) factor) will be the 
same for any complex ou. This behaviour near the ring singularity is universal, and thus 
independent of the un/stablc character of the field. 

To evaluate weather or not the above divergency is "reasonable" , we may compare with the 
static Maxwell field on Kerr that we get from the Kerr-Newman solution 

F^dA, A^ ^{dt- a siii^ (9) d(j)) (63) 

Note that since the Kerr-Newman metric is quadratic in Q, this field is a first order in Q 
solution of Maxwell equations on a fixed Kerr metric and, being Maxwell equations linear, 
(63) is also an exact solution on the Kerr background. For this static field, a straightforward 
calculation shows that ^ 

I static — 777 '■ 2)V4 (6^) 

2(r — lacosOp 

which exhibits the same degree of divergency as (62), the latter being possibly even milder 
along selected directions, or for some particular solutions. Note that the unstable solutions 
of the s = — 1 Teukolsky master equation, which evolve as e'^*/", decay in KNS as e"*^!^!/" for 
large |r|, as opposed to the slow, decay of the invariants of the static field abov. The 
unstable modes of Kill decay exponentially as r — >■ 00, and as a power of r — towards the 
inner horizon. 

In conclusion, we have shown that there are solutions of the Maxwell equations that behave 
in a similar way as the static field from the Kerr-Newman solution near the ring singularity, 
decay much faster away of along in spacelike directions, and grow exponentially with time. 



VI. CONCLUSIONS 

We have proved the existence of instabihties in both Kerr naked singularities (KNS) and 
the region beyond the inner horizon of sub extreme and extreme Kerr black holes (Kill). Our 
notion of instability is given by the existence of solutions of linear massless field equations 
that behave as e^*, A; > with a fast decay to zero as |r| — )■ 00 in the case of KNS (as 
r — )■ —00 and r — )■ r^^ in the case of Kill), where {t, r, 9, 4>} are Boyer-Lindquist coordinates. 
We have shown that there exist massless scalar fields, Weyl spinors. Maxwell fields and hnear 
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gravity perturbations with these properties. 

Since Kill and KNS do not admit a partial Cauchy surface, there is no natural notion of 
evolution from initial data, and therefore of instability in the usual way (bounded data grows 
unbounded). However, these spacetimes are time orientable, since the never vanishing vector 
field 

is everywhere timelike. The integral lines of V give a congruence of timelike curves filling 
the entire Kill (KNS) space which may be regarded as worldhnes of (accelerated) observers. 
Since the unstable fields we have found are axially symmetric (independent of 0), they 
grow boundless along these curves. That a congruence of observers measure a boundless 
growth of these linear fields (more properly, of any scalar made out from them) is to us an 
appropriate notion of instability in spacetimes lacking a partial Cauchy surface (for examples 
of evolution and stability notion in non globally hyperbolic spacetimes admitting a partial 
Cauchy surface sec the Schwarzschild and Reissner- Nordstrom cases [11], [12]).) 
A modification of the KNS metric based on alternative string motivated theories, or just 
the excision of a region around the ring singularity, gives rise to "superspinars" , compact 
rotating objects violating the black hole a < M bound [27]. The gravitational stability 
of these objects (more concretely, the region r > ro of KNS, assuming different boundary 
conditions at r = ro) was studied in detail in [28] and references therein. In Section 3.B 
of that paper, the superspinar gravitational instabilities results are confronted with those 
in our works [13] and [14], generalized in the present paper. This is done by studying the 
To — > — oo limit of unstable axially symmetric perturbations of a superspinar with a perfectly 
reflecting "string horizon" at tq — > —oo. This gives a result that agrees perfectly with our 
previous numerical calculations in [11]. 

The mathematical origin of the instability in both cases (KNS and superspinars) is, however, 
quite different. In the KNS case it is due to the negative portion of the effective potential 
(23) for small r < values, a region that is usually excised in a superspinar. The superspinar 
stability problem reduces also to a Schrodinger equation problem. However, in the KNS case, 
the Schrodinger equation is posed on the entire real axis, whereas in the superspinar case is 
a quantum mechanics problem on a half axis (which may well exclude the region where our 
effective potential is negative), with non trivial boundary conditions that are responsible of 
the instability. For Kill in the sub-extreme case we have also reduced the field equations to 
a QM problem on a half-axis, however, the boundary conditions are trivial, and the negative 
portion of the effective potential belongs to this half axis. 

We should also comment on the connection found in [28] (see also [29]) between circular 
geodesies with negative energies and superspinar instabilities. We note that those results do 
not seem to apply in a straightforward way to our case, since the arguments in [28] apply 
for £ = m >> 1 modes whereas unstable modes in this work have m = 0. We should also 
stress that, contrary to what is found in [28] the unstable modes whose existence is proved 
in this work exist for any value of a/M as long as a/M > 1. 
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Appendix A: Lower bounds to the radial potentials for the KNS as A; — > 0+ 

This Appendix gathers the calculations leading to the bounds given in (32) for the global 
minima of the potential (23) in the A; — )■ 0+ limit. These are used in Section III to obtain a 
bound if the radial spectrum in this limit. 



a. Case s = 



Introducing the adimensional variables x — r/M, a — a/M, the critical points of the 
potential (23) for s = are given by 



4r {a' + x^)(x^ - Zx' + a^x + a^) + a\a^ - 1)(1 - x) = 



(Al) 



For k small enough, three out of the five roots are real, and we order them as Xi < X2 < x^. 
Given the asymptotic behaviour (26), it is clear that the absolute minimum of V is reached at 
some of these points {xi or 2:3 if there are no inflection points). Inspection of the numerical 
solutions of (Al) for increasingly smaller k values suggests that Xi —00, X2 ^ 1 and 
X3 — )■ 00 as A; — )■ O"*". Guided by this observation we propose a solution of (Al) in the form of 
a power series in k taking the value X2 — 1 Bi, k — 0, and obtain by iteration the successive 
corrections in increasing powers of k, the result being. 



X2 



a\d' - 1) 



(A2) 



Then we assume x = x\\Ti (Al), solve for k. and expand the resulting expression for x\ ^ 00 
(which wc know that corresponds to taking k — )■ O"*") to extract the leading order behaviour of 
the relation between x\ and fc, which is x\ = — ((a^ — l)/4)^/'^/c~^/^. Inserting this expression 
plus a correction (to be obtained) back into (Al), we can iteratively get as many higher order 
terms as we wish. The first few of them are 



Xi = 



a 



1/4 



V2{Aa^ - 7) 
8(^2 _ 1)1/4 



k^/^ + 0{k). 



(A3) 



Proceeding in a similar way leads to 



X3 = 



V2{Aa^ - 7) 

8(q;2 _ 1)1/4 



k^/^ + 0{k) 



(A4) 



The value of the potential at these points are Vi — 1/4 + 0{k),V2 
1/4 + 0{k). It then follows that 



l/2 + C»(A;2),y3 



lim min{V {r,k,s = 0, a), r G M} = 1/4, 

fc-^0+ 



(A5) 



in agreement with (32) 
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b. Case \s\ = 1/2 



For s = 1/2, there are three real critical points xi < X2 < X3 {xi and 2:3 are local minima, 
X2 a maximum) only if k is bigger than an a dependent critical value, below which X2 and 
Xs coalesce into a single (inflection) point. In any case, the absolute minimum is reached at 
Xi, for which 

xi = -— - 1 ^—^ '-e + e» A;' , = k + + Oik"") 

2k a a 

For s = —1/2, the analysis is similar, with xi X2 and a global minimum at X3 — >■ 00 as 
k 0+: 

2A; q; cc^ a a"' 

It follows that 

lim min{y (r. A;, s = ±1/2, a), r e R} = 0, (A6) 

in agreement with (32) 



c. Case \s\ = 1 



For s = ±1, the critical points and potential values at these points are 



A(a^-3)k 8 (11 - 21) ^ 
6 a 9 



2^3 



the values of the potential at these points being 

Vi = -I ± -k + 0{k^) 
2 a 



V2^-^T -k + 0{k^) 
4 a 



Thus 



lim min{y(r, k,s = ±l,a),r eR} = -3/4, (A7) 

k-^0+ 



in agreement with (32). 
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d. Case \s\ = 2 

For s = ±2 the critical points and potential values are 

^^^^ 8(al-3)j 8 (47.^ 81)F ^ 
15 a 225 

A; 32 a 32 

1/3 

X3 



and and the values of V at these points are 

Vi = -l±-k + 0{k'); 
Z a 



It follows that 



4 a 
1 _ 315V3(c.^-l) 
^'"4 Ka2-l))V3^ +^^^ ^ 



lim min{1^(r, fc, s = 0, a), r G R} = -15/4. (A8) 



This completes the verification of (32). 



[1] Kerr, R. P., Phys. Rev. Lett 11, 237 (1963). 

[2] E. Newman and R. Penrose, J. Math. Phys. 3, 566 (1962); errata J. Math. Phys. 4, 998 (1963). 
[3] Teukolsky, S. A. Phys. Rev. Lett. 29 (1972) 1114; Astroph. J. 185, 635 (1973). 
[4] A. Komar, Phys. Rev. 113, 934 (1959). 

[5] O'Neill, Barrett (1995), The Geometry of Kerr Black Holes, Wellesley, MA: A. K. Peters. 

[6] J. L. Friedman, M. S. Morris, I. D. Novikov, F. Echcvcrria, G. Klinkhaninicr, K. S. Thorne, 
U. Ymtsever, Phys. Rev. D42 (1990) 1915-1930; S. W. Hawking, Phys. Rev. D46 (1992) 
603-611; J. L. Friedman and M. S. Morris, Coraraun. Math. Phys. 186, 495 (1997). 

[7] Penrose, R. in Black Holes and Relativistic Stars, ed. R. Wald, University of Chicago Press 
(1999); Earman, J. Bangs, Crunches, Whimpers, and Shrieks: Singularities and Acausalities 
in Relativistic Spacetimes, Oxford University Press (1995) and references therein. 

[8] Robinson, D.C., Phys. Rev. Lett 34, 905 (1975). 

[9] Mazur, P.O. in Proceedings of the 11th International Conference on General Relativity and 
Gravitation, ed. M.A.H. MacCahum, Cambridge University Press (1987) and references 
therein. [arXiv:hep-th/0101012vl.] 
[10] B. F. Whiting, J. Math. Phys. 30, 1301 (1989). 

[11] Dotti, G., Gleiser, R.J Class.Quant.Grav.23, 5063 (2006) [arXiv:gr-qc/0604021]; Class. Quant. 

Grav. 26 (2009) 215002. [arXiv:0809.3615 [gr-qc]]. 
[12] Dotti, G. and Gleiser, R.J, Class. Quant. Grav. 27, 185007 (2010) [arXiv:1001.0152 [gr-qc]]. 
[13] Dotti, G., Gleiser, R.J and Pullin, J. Phys.Lett. B644 (2007) 289-293 [arXiv:gr-qc/0607052]. 



22 



[14] Dotti, G., Gleiser, R.J, Ranea-Sandoval, I. F. and Vucctich, H. Class. Quant. Grav. 25, 245012 
(2008) [arXiv:0805.4306 [gr-qc]]; G. Dotti, R. J. Gleiser, J. Pullin, I. F. Ranea-Sandoval, 
H. Vucetich, Int. J. Mod. Phys. A24 (2009) 1578-1582. [arXiv:0810.0025 [gr-qc]]. 

[15] W. Kinnersley, J. Math. Phys. 10, 1195 (1969). 

[16] L. Kegeles, J. Cohen, Phys. Rev. D19 ,1641 (1979); P. Chrzanowski, Phys. Rev. Dll, 2042 

(1975). 

[17] E. Berti, V. Cardoso and M. Casals, Phys. Rev. D 73, 024013 (2006) [Erratum-ibid. D 73, 

109902 (2006)] [arXiv:gr-qc/0511111]. 
[18] Leaver, E. W. Proc. Roy. Soc. Lond. A 402, 285 (1985). 
[19] E. Seidel, Class. Quantum Grav. 6 (1989) 1057. 

[20] R. A. Brcucr, Gravitational Perturbation Theory and Synchrotron Radiation, Lecture Notes 

in Physics, Vol. 44 (Springer, Berlin, 1975). 
[21] R.A. Breuer, M.P. Ryan, Jr., and S. Waller, Proc. R. Soc. London A358, 71 (1977). 
[22] E. Berti, V. Cardoso and S. Yoshida, Phys. Rev. D 69, 124018 (2004) [arXiv:gr-qc/0401052]. 
[23] K. Meetz, II Nuovo Cimento XXXIV 690 (1964). 

[24] M. Reed and B. Simon, Methods of Modern Mathemaical Physics II: Fourier Analysis, Self- 

Adjointness, Academic Press, section X.l. 
[25] R. Wald, Phys. Rev. Lett. 41, 203 (1978). 

[26] O. J. C. Bias, H. S. Reall, J. E. Santos, JHEP 0908, 101 (2009). [arXiv:0906.2380 [hep-th]] 

[27] E. G. Gimon and P. Horava, Phys. Lett. B 672, 299 (2009) [arXiv:0706.2873 [hep-th]]. 

[28] P. Pani, E. Baraussc, E. Berti and V. Cardoso, Phys. Rev. D 82, 044009 (2010) 

[arXiv: 1006. 1863 [gr-qc]]. 
[29] Z. Stuchlik and J. Schee, Class. Quant. Grav. 27, 215017 (2010); S. R. Dolan and A. C. Ot- 

tewill. Class. Quant. Grav. 26, 225003 (2009) [arXiv:0908.0329 [gr-qc]]. 



23 



